The role of curvature in the slowing down acceleration scenario 
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We introduce the curvature 0.^ as a new free parameter in the Bayesian analysis using SNIa, BAO 
and CMB data, in a model with variable equation of state parameter w{z). We compare the results 
using both the Constitution and Union 2 data sets, and also study possible low redshift transitions 
in the deceleration parameter q{z). We found that, incorporating Q.k in the analysis, it is possible to 
make all the three observational probes consistent using both SNIa data sets. Our results support 
dark energy evolution at small redshift, and show that the tension between small and large redshift 
probes is ameliorated. However, although the tension decreases, it is still not possible to find a 
consensus set of parameters that fit all the three data set using the Chevalier- Polarski-Linder CPL 
parametrization. 

PACS numbers: 98.80. Cq 
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I. INTRODUCTION 

Since 1998 one of the biggest puzzles in cosmology is 
dark energy: the unknown component responsible for the 
current accelerated expansion of the universe 1] . In its 
simpler form, this can be described by a constant equa- 
tion of state parameter w = —I, corresponding to a cos- 
mological constant, leading to the ACDM model, the sim- 
plest model that fits a varied set of observational data. 

As the observations increase in precision, it is possible 
to reach levels that enable us to discriminate between a 
cosmological constant or a dynamical cosmological con- 
stant {w — w{z)). The source of this dynamical dark 
energy could be both, a new field component filling the 
universe, as a quintessence scalar field [2|, or it can be 
produced by modifying gravity l&\. 

In a recent paper [4| the authors suggest that cur- 
rent observational data could favor a scenario in which 
the acceleration of the expansion has past a maximum 
value and is now decelerating. The key point in deriving 
this conclusion is the use of the Chevalier-Polarski-Linder 
(CPL) parametrization [a], [Q| for a dynamical equation 
of state parameter 



w{a) — wq + {1 — a)wi, 



(1) 



where wq and wi are constants to be fixed by observations 
and a is the scale factor. An updated analysis performed 
using recent SNIa data was informed in [3] where similar 
conclusions were derived. Both analysis assumed a flat 
universe. 

Although it is often said that the inflationary scenario 
predicts a flat universe [3| , the situation has not been set- 
tled down yet, in part because the observations are not 
conclusive and because in theory flk = is not the only 
solution Iffl. In fact, the conclusion deduced from the 
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analysis of the observations using a w ^constant, imply- 
ing a model with cosmological constant |1 -I- wj < 0.06, 
only works if we assumed a flat universe. However, con- 
sidering the curvature as a new free parameter, it leads 
to an increment in the uncertainties in all the others pa- 
rameters. In particular, in the context of models with 
dynamical dark energy, a number of works have discussed 
the degeneracies between w{z) and flk [li| [I4] [13 UM 

It is the purpose of this paper to enhance the analy- 
sis of this phenomenon, relaxing the assumption of a flat 
universe, considering explicitly the geometric degeneracy 
between f2j. and w{z). This enable us to explore new fea- 
tures considering the CPL parametrization, enabling us 
to accommodate a similar behavior for q(z) using all the 
three observational test. This fact again shows the sen- 
sitivity of the observational data to even smaller values 
of the curvature parameter. The result found here shows 
that a single trend for q{z) is possible to draw using all 
the three observational probes, ameliorating the tension 
between data sets. 

The paper is organized as follows: in the next section 
we re-study the flat universe case using both SNIa data 
sets 16| , [Tg] . Then, we show how to generalize the anal- 
ysis considering a curved universe, in section III, which 
are our main results. We end that section with a discus- 
sion. 



II. A FLAT UNIVERSE ANALYSIS: UNION2 
VERSUS CONSTITUTION DATA 



The original analysis [J| was performed using the Con- 
sitution data set [I61 of SNIa. In this section we reanalyze 
this data set along the BAO and CMB constraints using 
the CPL parametrization. Then, we repeat the analysis 
using the full Union2 data sample [I3|, looking for dif- 
ferences in the consequences derived from these analysis. 
We end this section discussing our results along the ones 
obtained in the recent work Ml. 



The luminosity distance in a flat universe is 

dz' 



ddz) 



(1 + ^)^/ 



Eiz'V 



where 



E' 



'{z) = an(l + zf + (1 - r!,^)e'r ^^^^-'. 



(2) 



(3) 



The SNIa data give the distance modulus fiobs{z) as a 
function of redshift. We fit the SNIa with the cosmolog- 
ical model by minimizing the x^ value defined by 



XsNIa 



Np , , 



Mo6s(^i)]' 



(4) 
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where /i(z) = 5 log]^g[(iL (z)/Mpc] + 25 is the theoretical 
value of the distance modulus, and Np is the number of 
data points. 

The BAO data considered in our analysis is the dis- 
tance ratio obtained at z = 0.20 and z = 0.35 from the 
joint analysis of the 2dF Galaxy Redsihft Survey and 
SDSS data [l^, that can be expressed as 



Jv(0.35) 
Dv{Q.2Q) 



1.736 ±0.065, 



with 

Dv{zbao) 



zbao 
H{zbao) 



H{z) 



1/3 



(5) 



(6) 



We fit the cosmological model minimizing the x^ defined 
by 



Xbao 



[i:»y(0.35)/L»y(0.20) - 1.736]2 
0.0652 



(7) 



A result from the combination of SNIa and BAO is given 
by a joint analysis finding the best fit parameters that 

minimize X%Nia + X%AO- 

In addition, we incorporate to the analysis the CMB 
shift parameter [19|, which is the reduce distance at zis — 
1090 fH 



R 



yrt„iHir{zis) 



1.71 ±0.019. 



(8) 



Here we use the comoving distance r{z) from the observer 
to redshift z which is related to the luminosity distance 
by dL{z) = (1 + z)r{z). In this case we compute x^ by 



2 _ [i?-1.71]^ 

XCMB 



0.0192 



(9) 



to find out the result from CMB and the constraints from 
SNIa+BAO+CMB are given by XsNia + xIao + Xcmb- 
For the CPL parametrization, w — wq + wiz/{l + z), 
the second term in the righ hand side of ([3]) is 



(1 - a„)(l + z)3(i+-o+».) exp (-^) . ( 



10) 



TABLE I: The best fit values for the free parameters using 
the Constitution data set in the case of a flat universe model. 
See the related Fig. [T] 



Data Set 



Xmin 



n„ 



Wo 



Wi 



SN 461.23 0.4519 -0.2209 -11.228 

SN+BAO 461.60 0.4551 -0.1253 -12.255 

SN+BAO+CMB 466.94 0.2578 -0.9275 -0.019 



Using the Constitution sample consisting of 397 data 
points, assuming a fiat universe, leads to the results 
shown in Table HI They are very similar to the ones in- 
formed in [4] and [8]. Given the number of degrees of 
freedom N, a goo d est imate of the width of the x^ dis- 
tribution is cr = y^2/N, which in our case gives a ~ 0.07, 
a number taken into account to discriminate statistical 
significance among different runs. Having found these 
numbers, we can draw a profile of the deceleration pa- 
rameter q = —da/d^ as a function of the redshift z for 
each case. This is shown in FiglTJ Clearly, the incorpo- 
ration of the CMB constraints into the fitting process, 
leads to a completely different behavior for q{z). The 
values of Xmin show that the CPL parametrization is 
unable to fit the data simultaneously at low and large 
redshift. Because the best fit considering CMB gives a 
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FIG. 1: Using the Constitution data set [Ia| we plot the decel- 
eration parameter reconstructed using the best fit values for 
three cases: only SNIa (continuous line) , SNIa-|-BAO (dashed 
line) and SNIa+BAO-^CMB (dotted line). 

X^in 5 points higher than those considering SNIa and 
SNIa+BAO, indicates that CMB data favor a different 
behavior at small redshift, as the one observed in Fig IT] 
Considering the reduced x^ we notice that Xred = 1.169 
for SNIa+BAO and xled = 1-182 for SNIa+BAO+CMB, 
showing that using SNIa+BAO gives a much better fit to 
the model than the one using CMB, and again the inade- 
quacy of the CPL parametrization of fitting data at large 
and small redshift. The same can be concluded observing 
Figure [5] where the confidence contours in the wq — wi 
plane is plotted . Notice that the ellipses correspond- 
ing to SNIa alone and SNIa+BAO are almost identical, 
but the one considering CMB data does not overlap the 
others, indicating the tension between these data sets. 




q(z) 




FIG. 3: Using the Union 2 data set [IT]] we plot the decel- 
eration parameter reconstructed using the best fit values for 
three cases: only SNIa (continuous line), SNIa+BAO (dashed 
line) and SNIa+BAO+CMB (dotted line). 
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FIG. 2: Using the Constitution data set [la] we plot the la 
contours for the wq — wi parameters using the best fit val- 
ues for three cases: only SNIa (continuous line), SNIa-|-BAO 
(dashed line) and SNIa+BAO-fCMB (dotted line). 



Using the Union 2 set [l7|, the largest SNIa lumi- 
nosity distance sample currently available, consisting in 
557 type la supernovae, the analysis leads to the results 
shown in Table [III In this case, a ~ 0.06, then again 



TABLE II: The best fit values for the free parameters using 
the Union 2 data set in the case of a flat universe model. See 
also Fig. [1 



Data Set 



Xmin 



i 6ti 



Wo 



Wl 



SN 541.43 0.4197 -0.8632 -5.490 

SN-hBAO 542.11 0.4281 -0.7959 -6.537 

SN+BAO+CMB 543.91 0.2547 -0.9979 0.190 



the numbers quoted in the table have a trustable statis- 
tical meaning. However, by comparing this with the re- 
sults using the Constitution data set, we observe that the 
difference between the case with and without the CMB 
data is now smaller (~ 1.8), although still large com- 
pared to 3ct — 0.18. Plotting the deceleration parameter 
using these numbers leads to Figl^l Again, the case SNIa 
alone and SNIa-1-BAO are almost identical, but the case 
including CMB data does not shown the rapid change at 
small redshift as is also found using the Constitution data 
set. Similar results were recently published in i7|. Mean- 
while the statistical analysis using the Constitution data 
set give Xred values greater than one, indicating bad fit- 
tings, the ones obtained using the Union 2 sample, gives 
us over-fitting values, Xred ^ ^- This is important to no- 
tice, because an over-fitting work indicates that probably 
we are fitting noise instead of the actual relationship, and 
then more data would be necessary to improve the pre- 
dictive power of the model. Therefore, it is necessary to 
perform a study that compares both data sets of SNIa to 



shed some light on this point. 

This fact can also be viewed from Figure |4] where 
confidence contours in the wq — wi plane are plotted. 
Notice that again, the ellipses corresponding to SNIa 
alone and SNIa-|-BAO are similar, but the one consider- 
ing SNIa-l-BAO-l-CMB points into a completely different 
region. This confirm the previous claim that the CPL 
parametrization is not adequate here. 




Wo 



FIG. 4: Using the Union 2 data set [I3| we plot the la con- 
tours for the Wo — wi parameters using the best fit values 
for the three cases: only SNIa (continuous line), SNIa-|-BAO 
(dashed line) and SNIa+BAO+CMB (dotted hue). 



TABLE III: The best fit values for the free parameters using 
the Constitution data set in the case of a non flat universe 
model. See Fig. (5] 



Data Set 



Amin 



^ "Tl 



Wo 



Wl 



^k 



SN 459.15 0.7049 -0.4698 -3.6307 -0.9901 

SN+BAO 459.37 0.6861 -0.4868 -2.8019 -1.1534 

SN-f-BAO+CMB 461.14 0.4934 -0.1625 -11.2574 -0.0886 



III. RELAXING THE FLAT UNIVERSE 
ASSUMPTION 

The comoving distance from the observer to redshift z 
is given by 



r{z) 



r(z) 



c 1 



^ dz' 



sm 



{^/^n^r{z)) (11) 



Eiz'Y 



where ftk = —k/H^, k being the curvature constant and 



E'^iz) = nra{l + zf+nk{l + zf + nd,f[z), (12) 



f{z) = exp^3/ ^\l^'? dz' 



and Vide = 1 - ^m - ^fc- 

For the CPL paramctrization, w = wo + wiz / {1 + z), 
the function f{z) in p^ leads to 



/(z) = (l + z)3(i+-+-)exp(^-^ 



(13) 



Using the Constitution set we found the results displayed 
in table IIIII A direct comparison with the result of table 
|T1 where f2fc = was assumed, shows many interesting 
changes. In general, all the values for X^^^ ^'^s lower 
considering J7fc a free parameter. This result is expected, 
but the best fit values obtained are not. For example, 
in the cases where SNIa and SNIa-f BAO are considered, 
the best fit values are completely different: Vtk is around 
— 1, meanwhile the value for wi is four times less negative 
than in the flat case. Also the matter density becomes 
almost twice the value that in the flat case. 

Using the best fit values of the parameters we plot 
q{z) versus z for each case. We see from FigO that 
all the three cases behaves in a similar way. This is 
one of the main result of our work: the joint analysis 
SNIa-|-BAO+CMB using the curvature shows qualita- 
tively the same variation in q{z) for small redshift than 
those for SNIa alone or SNIa+BAO, a results that as far 
as we know has not been informed previously. 

The SNIa alone and SNIa-|-BAO look very similar to 
the already known, with the exception that q{z) is al- 
most reaching its zero value today z = Q. In each run, 
we have explored several different initial conditions, to 



q(z) 0.0 




FIG. 5: Using the Constitution data set [I3] we plot the decel- 
eration parameter reconstructed using the best fit values for 
three cases: only SNIa (continuous line), SNIa-(-BAO (dashed 
line) and SNIa-fBAO-fCMB (dotted line). 



TABLE IV: The best fit values for the free parameters using 
the Union 2 data set in the case of a non flat universe model. 
See Fig. [3 



Data Set 



/Cmin 



i Ln 



WO 



Wl 



Ofc 



SN 541.05 0.2973 -0.4362 -19.343 0.3551 

SN-fBAO 541.91 0.3433 -0.5488 -14.531 0.2499 

SN+BAO+CMB 542.18 0.4449 -0.8180 -5.1670 -0.0744 



find the best fit values. Regarding the case where is used 
SNIa-HBAO-HCMB, using a small value for Qk = -0.089, 
enable us to find a q{z) whose form follows the same trend 
that the previous cases. Here we must emphasize three 
important points: first, the present value of the decelera- 
tion parameter is positive today (z = 0), indicating that 
our universe is now decelerating, having crossed to pos- 
itive values at z ~ 0.1. Second, the form of q(z) shows 
that even considering the CMB data, the incorporation 
of the curvature to the analysis allows a transition in the 
deceleration parameter at small redshift z < 0.3. And fi- 
nally, although the SNIa-f BAO-fCMB case fit a q{z) that 
change sign, the profile of q{z) reaches a minimum value 
later (at z ~ 0.2) than in the other two cases (z ~ 0.4), 
crossing to negative values later than the other cases and 
also in the SNIa-|-BAO+CMB case the asymptotic value 
g ~ 1/2 for 2 > 1 is less than the one found in the other 
cases g ~ 1. However, this results does not imply that 
the tension between low and high redshift data has disap- 
peared. In fact, by drawing confidence contours around 
the best fit values in the wq — wi plane shown in figure 
ISl again the ellipses at one sigma does not overlap at all, 
even those for SNIa and SNIa+BAO. 

Using the larger Union 2 data set [17| we performed 
the same analysis as before obtaining the best fit param- 
eter shown in table IIVI The cases corresponding to SNIa 
data only and the joint analysis of SNIa and BAO data, 
behaves similarly to the ones already displayed using the 
Constitution data set as is evident in FiglT) However 
in this case the change in sign of q{z) take place later 
than the ones discussed in FiglS] Actually, this time 
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FIG. 6: Using the Constitution data set [la] we plot the la 
contours for the wq — wi parameters using the best fit val- 
ues for three cases: only SNla (continuous line), SNIa+BAO 
(dashed line) and SNIa+BAO+CMB (dotted line). 



q(z) 



U.b 




0.4 




0.7 










1 /' 


(HI 


1 A 




\ /' . 




\ /' • 


(17 


\ /' .' 




•A // . 


(1,4 






\ 11 ' 






0.6 





FIG. 7: Using the Union 2 data set [l3| we plot the decel- 
eration parameter reconstructed using the best fit values for 
three cases: only SNIa (continuos line), SNIa-|-BAO (dashed 
line) and SNIa+BAO-fCMB (dotted line). 



all the three fits behave similarly, reaching a minimum 
value at z ~ 0.2. Again, the value of the deceleration 
parameter to(iay is almost zero for the SNIa alone and 
SNIa-|-BAO cases, indicating that the acceleration of the 
expansion of the universe had already ceased, and in 
the SNIa+BAO+CMB case q is small but still negative. 
However, as was also the case using the Constitution set, 
the concordance in the profile of q{z) among the obser- 
vational probes is only apparent, because the one sigma 
confidence contours in the plane wq — wi indicate that 
although SNIa and SNIa+BAO share a region in that 
phase space, the CMB data spoiled the agreement. 

Comparing table HI with IIII[ and HIl and llVi we observe 
that the incorporation of the curvature parameter de- 
creases the value of Xmin ; indicating that a model with 
a non flat universe is preferred, even with a small value. 
Also is clear from the analysis that, both in the flat and 



non flat case, using the Constitution data set, leads to a 
wrong flt Xred > 1' meanwhile using the Union 2 data set, 
leads to a over- fitting, xled < 1' indicating the necessity 
of a further analysis of these SNIa data sets. 




FIG. 8: Using the Union2 data set |17'| we plot the la contours 
for the Wo — wi parameters using the best fit values for three 
cases: only SNIa (continuous line), SNIa+BAO (dashed line) 
and SNIa+BAO+CMB (dotted line). 



One of the main result from this analysis is that, the 
incorporation of the curvature as a new free parameter, 
enable us to accommodate a similar behavior for q(z) 
using all the three observational test. This fact again 
shows the sensitivity of the observational data to even 
smaller values (see Table |IV| of the curvature parameter. 
The result found here shows that a single trend for q(z) is 
possible to draw using all the three observational probes, 
ameliorating the tension between data sets. 

In this paper we have shown that incorporating the 
curvature fife as a free parameter in the Bayesian analy- 
sis using SNIa, BAO and CMB, enable us to ameliorate 
the tension between low and high redshift found in pre- 
vious works [J|,[3]- This is explicitly shown using both 
the Constitution [I^ and the recent Union 2 [131 data 
sets. Even considering small values for the curvature pa- 
rameter, it is possible to find a general behavior for q{z) 
at small redshift which shows a rapid variation in sign 
between z ~ 0.5 to z ~ 0. Because the CPL ansatz ([T]) 
implicitly assumes that nothing special occurs for DE in 
the redshift interval between z = to z ~ 2, the trend 
found in this work shows the incompatibility of the CPL 
parametrization in describing the variation of w{z) with 
redshift. This conclusion is reinforced based on the lack 
of overlap among the confidence contours in the three 
cases in the wq — wi plane. 
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